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Abstract: HEAT transfer by natural convection from a body to its finite enclosure is of importance 
in nuclear reactor technology, electronic instrumentation packaging, aircraft cabin design, the 
analysis of fluid suspension gyrocompasses, and numerous other practical situations. The steady 
natural convection heat transfer of fluids between two concentric isothermal spheres is investigated 
computationally with the help of FEV in ANSYS 14.5. The inner wall is subjected to a higher 
temperature and outer is at room temperature. The steady behavior of the flow field and its 
subsequent effect on the temperature distribution for different Rayleigh numbers and radius ratios 
are analyzed. 

Bossious boundary condition is taken for natural convection and which is solved in fluent 
module. Steady solutions of the entire flow field is obtained for Rayleigh number (5xl01<Ra<105), 
various Prandlt number and radius ratio (l.5<rr<3). The result shows that the Rayleigh number and 
radius ratio have a profound influence on the temperature and flow fields and Prandlt number has 
very negligible effect. The results of average Nusselt numbers are also compared with those of 
previous numerical investigations. Excellent agreement is obtained. 
Key Words: Convection, Nusselt Number, Ansys. 



HEAT transfer by natural convection from a body to its finite enclosure is of importance in nuclear 
reactor technology, electronic instrumentation packaging, aircraft cabin design, the analysis of fluid suspension 
gyrocompasses, and numerous other practical situations. Accurate prediction of such heat-transfer rates, which 
can now only be roughly approximated, is required in many engineering design problems. Since natural 
convection flow fields are buoyancy driven due to thermal effects, the thermal fields and hydrodynamic fields 
are very closely coupled, and knowledge of the flow field is essential to the complete understanding of the heat- 
transfer phenomena. Also, the unusual stability conditions associated with this problem are of fundamental 
interest in the fields of heat transfer and fluid mechanics. 

In 1968 analytical investigation of natural convection in annuli between two isothermal concentric 
spheres, the inner surface being hotter, was initiated by Mack and Hardee, the solution for steady axisymmetric 
natural convection between isothermal concentric spheres is obtained. 

In 1973 yin and powe conduct experimental investigation to describe concerning the natural 
convection flow patterns which occur in the annular space between two concentric isothermal spheres, the 
inner one being hotter. The several types of flow patterns observed are correlated with previously published 
temperature profiles and are categorized in terms of steady and unsteady regimes. 

In 1976 Douglass perform steady forced convection of a viscous fluid contained between two 
concentric spheres which are maintained at different temperatures and rotate about a common axis with 
different angular velocities is considered. The resulting flow pattern, temperature distribution, and heat- 
transfer characteristics are presented for the various cases considered. 

Ralph and scanlan 1977 results of a flow visualization study of natural convection in liquids contained 
between a heated sphere and its cooled cubical enclosure are reported Interference between up -flow and down- 
flow layers was noted for the smallest dimension ratio, and this created a rather unique flow pattern. 

Ramadhyani et al 1984 presents numerical finite difference solutions of combined natural and forced 
convective heat transfer in spherical annuli. The result the buoyancy effects can have a very significant impact 
on the heat transfer and fluid flow, particularly at low Reynolds numbers. 

Sanjay and sengupta 1988 examine the results of a numerical investigation of the natural convection 
process between isothermal vertically eccentric spheres with hotter inner core, and fount that Negative 
eccentricities have been found to enhance convection while positive eccentricities have the reverse effect. 
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I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 4 1 Iss. 101 Oct. 2014 1291 



Computational Analysis of Natural Convection in Spherical Annulus Using FEV 



Results also show that heat transfer actually increases slightly for very high positive eccentricities where 
conduction plays an important role. 

Garg in 1992 give a finite-difference solution for steady natural convective flow in a concentric 
spherical annulus with isothermal walls has been obtained. The stream function -vorticity formulation of the 
equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. 

Chiu and Shich 1999 perform A numerical analysis determine the heat transfer of micropolar fluids 
by natural convection between concentric spheres with isothermal boundary conditions Results indicate that the 
heat transfer rate of a micropolar fluid is smaller than that of a Newtonian fluid, and the main controlling 
parameter is the dimensionless vortex viscosity. 

Vadim and rath 2002 The energy stability problem with respect to axisymmetric disturbances of the 
natural convection in the narrow gap between two spherical shells under the earth gravity is discussed. The 
results are compared with the results of the linear stability analysis for the same problem. 

Wen and Yen 2008 the effects of micro-rotation and vortex viscosity in micropolar fluids have been 
investigated numerically to determine heat transfer by natural convection between concentric and vertically 
eccentric spheres with specified mixed boundary conditions. 

InChen 2010 investigates the effects of eccentricity and geometric configuration with a Newtonian 
fluid, numerically to determine heat transfer by natural convection between the sphere and vertical cylinder 
with isothermal boundary conditions. Results of the parametric study conducted further reveal that the heat and 
flow fields are primarily dependent on the Rayleigh number, eccentricity and geometric configuration, for a 
Prandtl number of 0.7, with the Rayleigh number ranging from 103 to 106, the three eccentricities and two 
geometric configurations. 

Alassar 2011 An exact solution of the problem of heat conduction in the annulus between eccentric 
isothermal spheres with internal heat generation is obtained. The solution is given in terms of the temperature 
distribution and local and average Nusselt numbers. 

Sangita 2013 reports results of a numerical investigation of natural convection in a spherical porous 
annulus. The inner and outer surfaces are subjected to constant temperatures. The Brinkman extended Darcy 
flowmodel is considered in her study. 

Hatami and ahangar 2014 present, temperature distribution equation for a fully wet semi -spherical 
porous fin. The driving forces for the heat and mass transfer are considered temperature and humidity ratio 
differences, respectively. , the effects of porosity, Darcy number, Rayleigh number, Lewis number, etc. on the 
fin efficiency are investigated. 

II. Mathematical Modeling 




Figure 1 Physical model and coordinate system and Sectional view 



figure 1 Here 9 is the polar angle and varies from 0 to n The governing equations in polar coordinate 
are described in dimensional form specifying boundary conditions for velocity and temperature as below: 

— (r 2 u sin 0) + — (r 2 v sin 0) = 0 
dr 36 
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Boundary Condition 

u=0, v=0, T=Ti at r=ri for 0<6 
u=0, v=0, T=Ti at r=ro for 0<6<7T 
u=0, v=0, dt/d9=0 at 9=0 for ri<r< ro 
u=0, v=0, dt/d9=0 at 9=7i for ri<r< ro 

Non Dimensional Parameters 
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Figure 2 Model Geometry 
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Figure 3 Mesh Model 



III. Results 

Natural convection between two concentric spheres, inner at constant higher temperature while the 
outer is at low temperature, filled with a fluid has been studied. Now the various results are obtained from the 
analysis using FEV scheme. The various effects of three non-dimensional parameters Rayleigh number (Ra), 
Prandlt number (Pr) .and radius ratio(rr) on flow field, temperature field and Nusselt number have been 
discussed in this chapter. 

A. Grid Sensitivity Test 

Before studying the influence of the above three parameters on temperature field and flow field, the 
optimum values for some parameters used during calculation by FEV. Those parameters are optimum Grid 
size, acceleration factor (co) and the optimum error tolerance limit (e). First to the find the optimum grid size 
for numerical result the grid sensitivity test was conducted as shown in the table 1 



Table 1 Table Sensitivity test for rr = 2, Pr = 0.7 and Ra = 10 3 



Grid Size 


Nu h 


Nu c 


Avg. Nusset Number 


20x20 


1.99 


1.99 


1.99 


30x30 


1.97 


1.97 


1.97 


40x40 


1.96 


1.96 


1.96 


50x50 


1.96 


1.96 


1.96 



Examining Table- 1 shows that the smaller grid size 40x40 and 50x50 did not cause any significant 
changes in the magnitude of the average Nusselt number for the test cases a uniform grid size of 40x40 was 
chosen for calculation of all cases in the study. 

On the basis of the grid size 40x40 the optimum value of acceleration factor was checked According 
to SAR scheme in literature the acceleration factor co varies from 0 to 2. 

Acceleration factor co = 1 has been used for all the calculations having Ra < 3x 10 4 while acceleration 
factor co varies from 0.5 to 0.1 for higher Rayleigh due to the convergence criterion. 

Table 2 shows the values of average Nusselt number for different error tolerance limit. Nusselt 
number value increases as the value of error tolerance limit increases and attains a steady state at s=10" 5 . As 
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there has been marked very less changes in value between E value 10-5 and 10-6 in order to reduce the CPU 
time 10" 5 has been used as optimum error tolerance limit. 



Table 2 Optimum error tolerance limit (e) for Ra= 1000, Pr= 0.7 an d rr= 2 



Eps. Value 


Nuh 


Nuc 


Avg.Nu 


0.0001 


3.40459 


3.52778 


3.466185 


0.00001 


.4281 


3.45719 


3.442645 


0.000001 


3.42988 


3.45429 


3.442085 



All these values in table 2 were calculated for grid size 40x40, Ra= 10 5 , Pr=0.7, rr=2, co t = co 0 = co s 
=0.5. While for table -3 Ra= 10 4 , Pr= 0.7, rr=2, co t = co 0 = co s =1. 



B. Validation of FEV scheme 

The values of different parameters like average Nusselt number and maximum stream function 
calculated using FEV and SAR scheme were being validated by earlier papers as seen in Table 3. 



T able 3 Comparison of Present work with reference papers at rr= 2, Pr= 0.7, Ra= 10 3 . 



Reference 


Avg. Nu 


max 


H.W. Wu, Wen, C.Tsai [10] 


1.10 


3.25 


H.S. Chu, T.S.Lee [7] 


1.09 


3.21 


Mack and Harde [1] 


1.12 


3.21 


Present 


1.1 


3.23 



Table.4 compares the results of present work with results from different paper, which shows the 
results are very accurate. The percentage of error varies between 0.90-1.181 for average Nusselt number and 
between 0.62- 7.38 for maximum stream function. 

In Table 4 comparison of the present values has been done with the values as given in the paper of 
H.S. Chu and T.S.Lee[8] for different Rayleigh number. Before calculating the present data's for comparison, a 
relation is being developed between the expression for Rayleigh No used by Chu and Lee and the Present 
Rayleigh No. As it can be noticed from the table. 5 that all our values are very much accurate. The relation 
being developed is 



Ra p = (rr - 



Where, Ra = 



Ra p = 



if x Ra 
3a 

gj3AT(r o - ri ) 3 



expression of Rayleigh number used by H.S. Chu and T.S.Lee 



9a 



expression of Rayleigh number used in present work 



Table-4 Comparison of FEV and SAR scheme result with H.S.Chu and T.S.Lee at Pr= 0.7 for different 

Ra value and radius ratio 



rr 


Ra 


Ra p 


Nusselt number 
(H.S.Chu and 
T.S.Lee) Ref. 
[10] 


Average 
Nusselt 
number 
Present 


1.2 


lxlO 3 


8 


1.00 


1.00 


1.2 


lxlO 4 


8X10 1 


1.00 


1.00 


1.2 


lxlO 5 


8xl0 2 


1.01 


1.01 


1.5 


lxlO 2 


1.25X10 1 


1.00 


1.00 


1.5 


lxlO 3 


1.25xl0 2 


1.00 


1.00 


1.5 


lxlO 4 


1.25xl0 3 


1.07 


1.07 


1.5 


lxlO 5 


1.25xl0 4 


1.92 


1.95 


1.5 


lxlO 6 


1.25xl0 5 


3.71 


3.47 


2.0 


lxlO 2 


lxlO 2 


1.00 


1.00 


2.0 


lxlO 3 


lxlO 3 


1.10 


1.10 


2.0 


lxlO 4 


lxlO 4 


1.97 


1.96 


2.0 


lxlO 5 


lxlO 5 


3.49 


3.43 
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C. Flow and Temperature fields 

Fig. 4 and Fig. 5 present streamline configurations, isotherms for radius ratio of 2.0, Prandlt number 
0.7 at two different Rayleigh number 10 4 and 10 5 . These results are designed to show the influence of the 
Rayleigh number on the flow field and isotherms. The fluid in the close vicinity of inner sphere has a lower 
density than that of near the outer sphere because the inner sphere is at higher temperature than that of outer 
sphere. Thus the fluid near the surface of the inner sphere moves downward. As the fluid moves downward, it 
loses energy and eventually forces the separation of the thermal boundary layer along the outer sphere. The 
heavy fluid then enters the thermal boundary of the inner sphere and completes the recirculation pattern. The 
centre of crescent shaped eddy stayed close to midgap but moved into the upper hemisphere as the Rayleigh 
number increases. This is apparent as we compare the radius ratio 2.2 in the Fig.4, where the Rayleigh number 
is 10 4 , with that in the Fig. 5, where Ra=10 5 . For the case of Ra= 10 5 laminar convection is the more dominant 
of heat transfer than in Ra=10 4 As the circulation of fluid is more dominant in the upper zone making the 
outer layer warmer, while the lower domain is stagnant. The transport of the hot fluid to the outer layer is more 
prominent from the Fig. 4 and 5. For the conditions selected for fig.4 and 5 Pr=0.7 at different Ra value at rr 
value 2.2 the maximum values for stream function moves up ward as the Ra value increases. This behaviour 
can be shown from Fig 4 where for Ra=10 4 \|/ max = 14.138164 obtains at angle =63°, while for Ra=10 5 28.772739 
at angle =58.500004°. 




Figure 5 Isotherms and stream lines for Pr=0.7, rr = 2.2 at Ra= 10 5 for Velocity distribution 
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Figure 6 Isotherms and stream lines for Pr=0.7, rr = 2.2 at Ra= 10 4 for Static temperature 




Figure 7 Counter of Stream Function of Spherical Annulus 




Figure 8 Isotherm and stream line contour of Static Pressure 
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Contours ofWall Shear Stress (pascal) 

Figure 9 Isotherm and stream line contour of Wall shear stress 
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Figure 10 Variation of Nusselt number at different Pr values for Ra= (10 - 5x10) at radius ratio 2. 
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Figure 1 1 Variation of average Nu number and Ra. Number (5x10 -10 ) at different Prandlt No. for radius 

ratio 2 
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Figure 10 and 11 shows the variation of avg. Nusselt Numbers with Ra number at different Prandlt 
No. for radius ratio (rr = 2). Fig. 10 shows the variations of avg. Nu. No where Ra. No. varies from lO^x 
10 3 .and in Fig 11 Ra. No. varies from 5xl0 3 -10 5 . From the both the graphs two characteristics are being 
noticed; there is no effect or slight effect of Pr. No. on the heat transfer as the Rayleigh No increase. Another 
one, the Nusselt No varies almost linearly with Rayleigh No. As in the graphs the lines are nearly straight. 



.*.5 




1.0 1.5 2.0 2.5 3.0 3.5 

Radius Ratio, rr 

Figure 12 Variation of average Nusseltnumber and Radius Ratio (rr) at different Rayeigh value. 

Figure 12 shows the variation of avg. Nusselt No and the Radius Ratio at different Rayleigh No. Here 
we have taken Pr= 0.7 (for air) fixed. It has been seen that there is no influence of Prandlt values up on Nusselt 
number variation. It can also be conclude that for every Ra. value there is a critical radius ratio. Up to which 
the value of Nusselt No increases and after that it starts decreasing. The value of critical radius ratio decreases 
as the value of Rayleigh No. increases. Both the value of critical rr increases when Rayleigh No. varies up to 
103. 
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Figure 13 Steady local Nusselt number versus angular position at radius ratio 2 and Pr=0.7 at different Ra 

number 

Fig. 13 shows the distributions of steady local Nusselt numbers related to angular position when 
Rayleigh number changes, the effect of different Rayleigh number can be seen. According to the conditions 
Pr=0.7, rr =2, the local Nusselt number at the exterior of inner sphere have minimum value at a~oO.This 
implies that at upper symmetry the inner sphere has conduction and as the angular position increases natural 
convection becomes more dominant. While on the contrary at 6=0° the portion near the outer sphere has 
convection as dominant mode for heat transfer and as the angular value increases the mode of heat transfer 
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changes towards conduction. The value of Nuh increases, the maximums occurs at 9=180°. The local Nusselt 
numbers on interior of outer sphere (Nuc) decreases gradually as angular position increases. The minimums 
are occurred at 9=180°. The whole variation is larger than Nuh. Nuc, decreases because the heat boundary 
layer becomes thicker as it moves downward and generates stagnation area at the bottom of outer sphere. 
Therefore the Nu, have big variation with angular positions as shown in fig. 13. 

IV. Conclusion 

The present work investigates computational the natural convection heat transfer of an Newtonian 
fluid contained between two concentric isothermal spheres. The behaviour of flow field and temperature field 
are shown graphically. The present results of local and average Nusselt number are in good agreement with the 
earlier works. The major results may be summarized below as follows: 

At fixed radius ratio the average Nusselt number increase with increase in Rayleigh number. 

At fixed Rayleigh number the Nusselt number first increases with increase in radius ratio, but after a 

critical radius ratio value Nusselt number stars decreasing. 

The centre of main vortex, where stream function is maximum, moves upward with increase in both 
Rayleigh number and radius ratio. 

Prandlt number has very slight effect upon Nusselt number if radius ratio has kept constant. 
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Nomenclature 



g 


acceleration due to gravity thermal diffusivity of the fluid, [m/s 2 ] 


k 


Thermal conductivity of the fluid the fluid, [W/m-k] 


Pr 


Prandtl number of the fluid 




Radii of inner and outer spheres respectively, [m] 


d 


thickness, [m] 


rr 


Radius Ratio 


R 


Dimensionless radial coordinate 


e 


Dimensionless angular coordinate 


Ra 


Rayleigh number based 


Gr 


Grashof number 


T 


Temperature of the fluid [K] 


Ti, To 


Temperature of the inner and outer spheres, respectively[K] 


u 


velocity component in the r- direction, [rn/s 2 ] 


U 


dimensionless velocity component r- direction 


V 


velocity component in 8 -direction, [m/s 2 ] 


V 


dimensionless velocity component in in 8 -direction 


Nu h 


Average Nusselt Number on hot sphere 


NU C 


Average Nusselt number on cold sphere 


AR 


Small interval in R-direction 


A9 


Small interval in B direction 



Greek Symbols 



a 


Thermal diffusivity of the fluid 


P 


Coefficient of volumetric expansion 




Kinematic viscosity 


P 


Density of fluid 


pm 


Density of fluid at Tm 


AR 


Mesh size in the radial direction 


A9 


Mesh size in the B -direction 


T 


Dimensionless temperature 


Y 


Kinematic viscosity of the fluid 




Stream function 




Vorticity 
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